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Abstract. A theory is presented for the behavior of an array of multi-lamellar vesicles (the onion phase) 
upon addition of solvent. A unique feature of this system is the possibility to sustain pressure gradients by 
tension in the lamellae. Tension enables the onions to remain stable beyond the unbinding point of a flat 
lamellar stack. The model accounts for various concentration profiles and interfaces developing in the onion 
as it swells. In particular, densely packed 'onion cores' are shown to appear, as observed in experiments. 
The formation of interfaces and onion cores may represent an unusual example of stabilization of curved 
interfaces in confined geometry. 

PACS. 83.70.Hq Heterogeneous liquids: suspensions, dispersions, emulsions, pastes, slurries, foams, block 
copolymers, etc. - 87.16.Dg Membranes, bilayers, and vesicles - 82. 65. Dp Thermodynamics of surfaces 
and interfaces 



1 Introduction 



The self-assembly of amphiphilic molecules (surfactants) 
in solution has been the subject of intensive research in the 
past decades fjl@]. Such molecules may self-assemble into 
a wide variety of aggregate morphologies such as micelles, 
bilayers, bilayer stacks (the lamellar, or L a phase) and 
various other liquid-crystalline structures. The richness of 
these phenomena is both useful for numerous chemical 
and biochemical applications and challenging for basic re- 
search. 

Yet another morphology of amphiphilic self-assembly 
has recently drawn considerable attention. When a lamel- 
lar (L a ) phase is subjected to shear, it undergoes a dy- 
namic transition into an array of close-packed multilayer 
vesicles, referred to as the onion phase y,Q.g|jy,|3l • The typ- 
ical size of an 'onion' is of order 1-10 /im, whereas the 
inter-membrane spacing is of order 10 nm, and the thick- 
ness of the membrane itself is of order of a few nm. Thus, 
each onion comprises a spherical stack of more than sev- 
eral hundred spaciously packed membranes. Although this 
structure is not an equilibrium one, it may remain stable 
for days. In addition, one can control the size of the onions 
by changing the shear rate || , and encapsulate small par- 
ticles inside them g. These appealing features have po- 
tential applications, e.g., in the pharmaceutical industry. 
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Various properties of the onion phase are not well un- 
derstood theoretically. In particular, the detailed mecha- 
nism of the L^-to-onion dynamic transition has not yet 
been established M. Simple scaling arguments have been 
successful in accounting for the dependence of the onion 
size on shear rate and inter- membrane spacing [p|,|l0f , as 
well as the viscoelastic behavior of the onion phase |11| . 

In the current work we focus on another property of 
the onion phase — its behavior under dilution. This aspect 
has been investigated in recent experiments |l2Ulq] ■ When 
a regular L Q phase is diluted, the added solvent is accom- 
modated in the inter-membrane spacings. The membrane 
stack thereby swells, until a 'melting' transition into a dis- 
ordered L 3 ('sponge') phase occurs 0,H- The onion phase 
is observed to dissolve into L3 as well. Yet, because of the 
confined spherical geometry, the stack cannot swell with- 
out a supply of additional surfactant. (Amphiphilic mem- 
branes, to a good approximation, are practically inextensi- 
ble |[I5[ .) As a result, the dissolution progresses very slowly 
(over days) through onion coalescence and disintegration, 
giving rise at intermediate times to various defects and 
instabilities. In particular, the formation of dense onion 
cores has been observed in two different systems p2| , [l3| , 
as is demonstrated in Fig. pi 

Our model of onion swelling relies on an assumption 
concerning separation of time scales. The time scale of 
swelling of the entire onion is assumed to be much longer 
than the one required for internal equilibration within the 
onion. Experimentally, the former is of order of a few 
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hours. The latter is related to the time scale character- 
izing the formation of passages ('necks') between adjacent 
membranes flq] . We thus assume that passages form on 
time scales much shorter than hours. Though reasonable, 
the validity of this assumption is still to be established ex- 
perimentally. There is a third time scale, corresponding to 
onion coalescence and disintegration into L3. Experimen- 
tally, it is found to be of order of days and will be ignored 
here. Thus, although the onion phase is a system far from 
equilibrium, the above assumption allows us to present an 
essentially equilibrium model of onion swelling. 

A fundamental question arising from the study of onion 
dissolution concerns the possibility to stabilize a curved 
interface between two coexisting domains in confined ge- 
ometry. In common situations such a possibility does not 
exist. (For example, a liquid droplet inside a coexisting 
vapor phase is unstable, and will either shrink and vanish 
or expand to a macroscopic phase |17| ; this is the origin 
of nucleation barriers in lst-ordcr phase transitions.) Sta- 
ble interfaces in confined geometry are found in specific 
systems, such as magnetic garnet films, block-copolymer 
melts, phospholipid monolayers and microemulsions, ex- 
hibiting modulated phases pq| . Competing interactions in 
such systems lead to a negative surface tension (i.e., neg- 
ative stiffness term in a coarse-grained Ginzburg-Landau 
free energy), which is stabilized by higher-order terms. 
This gives rise to a finite characteristic length scale of in- 
terface modulation p8| . As seen in Fig. [j], onions under 
dilution seem to exhibit a stable interface between a con- 
fined, dense core and a dilute shell. We shall try to demon- 
strate that this behavior might represent a new way to sta- 
bilize a curved interface, arising from the unique ability of 
onions to sustain pressure gradients at equilibrium. 

In section |2| we extend a simple theory for the unbind- 
ing transition in a flat L Q phase [ |19| to the case of mem- 
branes with tension. Based on this extension, we formulate 
in section H a model of a single onion as an inhomogeneous, 
spherical membrane stack. We then present the resulting 
concentration and tension profiles in the onion. Finally, in 
section ^, we summarize the results and point at future 
directions. 



2 Lamellar Phase with Tension 

Fluid membranes in a lamellar stack experience steric re- 
pulsion arising from their reduced undulation entropy (the 
Helfrich interaction) |20) . In the case of tensionless mem- 
branes the interaction energy per unit area is 



fund(D) 



bT 2 

kD 2 ' 



(1) 



where T is the temperature (in energy units, i.e., fee = 1), 
k the bending rigidity of the membranes, and D the inter- 
membrane spacing. (The numerical prefactor b is still un- 
der controversy plj |; Helfrich's calculation pOfgives b = 
37r 2 /128 ~ 0.2, whereas computer simulations |22| yield a 
lower value of b ~ 0.06.) The long range of the Helfrich in- 
teraction, Eq. (|lj), results from the 'Soppiness', i.e., strong 



thermal undulations, of tensionless membranes. The inter- 
play between the Helfrich repulsion and other, direct in- 
teractions determines when a lamellar stack of membranes 
becomes unstable and unbinds. The rich critical behavior 
exhibited by this system was thoroughly investigated us- 
ing functional renormalization group techniques 1 23 04J . 

Subsequently, a much simpler theory for the unbind- 
ing of a lamellar membrane stack was proposed ||19| . It 
employs a similar argument to Flory's for polymers - 
since the stack is a 'soft', entropy-dominated system, one 
has to accurately account for entropy [i.e., the Helfrich 
repulsion, Eq. (Q)], while the other interactions can be 
incorporated in an approximate, 2nd-virial term. The re- 
sulting (grand-canonical) free energy per unit volume of 
the lamellar stack is 



1 



f(<f>) = 7iF - XF - M 



(2) 



where <f> = 5/D is the surfactant volume fraction (<5 being 
the membrane thickness), \ is a 2nd-virial coefficient, and 
fj, the surfactant chemical potential. All energy densities 
have been scaled by 2bT 2 /(n5 3 ). For \ > the free en- 
ergy ([2|) describes a Ist-order unbinding transition as /x 
is lowered. The chemical potentials and volume fractions 
corresponding to the binodal and spinodal of this transi- 
tion are 
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The free energy (g) also has a critical point at \x = x = 
0, which is of much theoretical interest p3[ , but of no 
relevance to the current discussion; hereafter, a positive 
value of x is assumed. 

Consider now a stack of membranes having tension 
a. Surface tension strongly suppresses membrane undula- 
tions and, hence, the fluctuation-induced interaction be- 
tween tense membranes has a much shorter range. The 
calculation of this interaction is more complicated than for 
tensionless membranes plfl . Renormalization-group calcu- 
lations p5| and computer simulations |26[ yield an expo- 
nential decay with distance. A simpler, self-consistent cal- 
culation of this interaction P7J gives 



/u„d(£>) 



bT 2 
kD 2 



D/h 



smh(D/l T ) 
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where It is the length arising from the combination of 
tension and thermal energy, It = (2T/na) 1 / 2 . The di- 
mensionless parameter x = D/It — S/(It4>), depending 
on both a and </>, determines whether the tension has a 
significant effect on the interaction. For 1 C 1 fund co- 
incides with the tensionless expression, Eq. (|lj), whereas 
for x ^>> 1 the tension strongly suppresses membrane un- 
dulations and the interaction decays exponentially with 
distance, in accord with the renormalization-group result 
[ |25| . (The numerical prefactor in It was chosen so as to 
recover the renormalization-group result for high tension.) 
For brevity we use hereafter the following notation: 



G(x) 



1 



2 sinh x 



, . dG cosh a; 

9( x ) = ~r = . ; 3 ■ 
ax sinh x 
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, dg 4 sinh x — 3 

ax smh a; 

The 'Flory-like' free energy per unit volume in the 
tense case is 



f{ ( t>,a)=- X 2 G{x)f-xct> 2 -^ 



(5) 



where the energy densities have been scaled, again, by 
2bT 2 /(nS 3 ). (Note that x should not be taken as an inde- 
pendent degree of freedom but as (^-dependent. The two 
independent degrees of freedom are (f> and a. Nevertheless, 
x is extensively used below, in order to make the formu- 
lation more concise.) 

Due to the additional degree of freedom — membrane 
tension — bound stacks can be stabilized even beyond the 
unbinding point of the tensionless case, i.e., for \i < [i sp . 
Hence, instead of a transition point there is a transition 
line, a* (/i), such that the stack is bound for a > a* , and 
unbound for a < a* . The equations for the binodal and 
spinodal lines are: 



binodal: 2 
spinodal: 6 



G(x) + xg(x) [i 

[x 2 g(x)] 2 jUbin ' 

G(x) - xg(x) - x 2 g'(x) 



X 3 g(x) 



g'{x)f 



JL 

Msp 



2,\ 



x 4 g'(x) 



(6) 



The two lines are drawn in Fig. J2|A. Figure |p shows the 
line x*(/j.) corresponding to the binodal. When \i is slightly 
lower than //bin, the value of x* required to stabilize the 
stack increases sharply. On the other hand, values of x* 
much larger than 1 are required only for chemical poten- 
tials much lower than /ibin- Hence, realistic values for x* 
should be of order 0.1-3. 



3 Onion Swelling 

3.1 The Model 

In view of the assumption regarding separation of time 
scales, presented in section [l], we consider the onion as a 
spherical stack of constrained size. (Recall that the initial 
size of an onion is, by itself, determined not by thermo- 
dynamic equilibrium but by the shear rate that led to 
its formation J5|.) As dilution progresses, the constrained 
radius, R, increases slowly, such that at any instant the 
interior of the onion can be assumed in thermodynamic 
equilibrium. This implies, in particular, that the entire 
onion has a single, uniform chemical potential, \x. Since 
the temperature, T, and the total number of surfactant 
molecules in the onion, N, are taken as fixed, \i must con- 
tinuously decrease upon swelling. Thus, the equilibrium 
ensemble relevant to the actual dilution process (at time 
scales shorter than hours) is that of fixed (T, R, N) , where 
R is regarded as a slowly increasing external parameter, 



leading to a slow decrease in \x and po, the osmotic pres- 
sure. For mathematical convenience, the model is formu- 
lated in the equivalent ensemble of fixed (T,po,/j,), where 
fj, is is a slowly decreasing parameter. 

As mentioned in section FP, each onion comprises a 
stack of many spaciously packed membranes. It is there- 
fore justified to employ a coarse-grained, continuum model. 
Unlike a flat lamellar phase, the inter-membrane spacing 
in the onion is not expected to be uniform. Hence, we al- 
low for non-uniform profiles, writing a Ginzburg-Landau 
free energy of the form 



F[4>{T),a(r)]= /dr 



f?|V0| 2 + /(0,a) 



-in 



-PoR" 



The free energy density, f(<f>, a), has been defined in equa- 
tion (ph, and the integration is over the volume of the 
onion. An energetic penalty for concentration gradients 
has been included in equation fm), where D is a stiffness 
coefficient. The Lagrange multiplier po should ensure that 
the onion radius has the constrained value R. 

The state of the onion is defined by the concentration 
profile 4>{r) (or, equivalently, by the set of inter-membrane 
spacings). The equilibrium concentration profile is thereby 
found from a variation principle. A major theoretical com- 
plication is the fact that the tension in the membranes 
enters at two distinct levels — as a macroscopic parameter, 
e.g., balancing the pressure difference across a membrane, 
and as a microscopic parameter having a dramatic effect 
on membrane interaction. In order to overcome this obsta- 
cle we use a self-consistent scheme. We assume that the 
membranes have attained certain values of tension, lead- 
ing to a (still unknown) profile x(r). Subsequently, we find 
the resulting concentration and pressure profiles and then 
require self-consistency, i.e., that the presumed tension in 
the membranes balance the pressure gradients. Note that 
this self-consistency does not require that the tension pro- 
file be smooth. Hence, no penalty for spatial changes in 
tension has been included in equation M); near-by mem- 
branes may have very different tensions. In practice, the 
formation of passages should act to equalize the tension 
between layers and, hence, a certain penalty for tension 
differences is expected. We assume that such a tension- 
gradient term would not have a drastic effect on the re- 
sults. Hence, we omit it to avoid addition of another pa- 
rameter to the model. 



3.2 Profile Equations 

By taking the variation of F with respect to <f>{r < R) we 
obtain the first profile equation, 

f2V 2 (j) + x 2 [G(x) - xg{x)](j) 2 +2 X (j) + H = 0, (8) 

where, for our spherical-symmetric case, V 2 = d 2 /dr 2 + 
(2/r)d/dr. (Recall that the variation is taken while keep- 
ing <7, not x, fixed.) 

Equation (@), obtained from a variation of a Ginzburg- 
Landau functional, is of the generic form widely used to 
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study interfaces p8[ . It is equivalent to imposing a uniform 
chemical potential throughout the system. Such a profile 
equation is usually supplemented by boundary conditions 
for the order parameter and its gradient far away from the 
interface, that ensure the uniformity of pressure. Such a 
system of equations is over-determined (a 2nd-order equa- 
tion with four boundary conditions). It has a solution for 
a flat geometry, but does not have one for a confined (e.g., 
spherical) geometry. Thus, such a Ginzburg-Landau for- 
malism cannot in general produce stable, confined inter- 
faces. However, the onion phase has a very special prop- 
erty. Since the system is composed of concentric closed 
sheets, the pressure need not be uniform throughout the 
system; pressure gradients can be sustained by the tension 
in the membranes. In particular, if the system is divided 
into coexisting domains, the pressure does not have to 
be equal in the different domains. This property bypasses 
the usual uniform-pressure boundary conditions. They are 
to be replaced by an equation, derived below, balancing 
pressure gradients and tension. 

Using Green's identity and equation (g), we can rewrite 
the free energy at equilibrium as 



where e = n/^TrbT). 

The profile equations (||) and (|l2|), along with appro- 
priate boundary conditions, determine the profiles <p(r) 
and x(r). 



3.3 Boundary Conditions 

The onion is in contact with the surrounding environ- 
ment through its outer layer. This layer requires a sepa- 
rate treatment, so as to yield the boundary conditions for 
the profile equations (g) and (|l2|), i.e., 4>(R), x(R) and 
d(f>(R)/dr. We employ a self-consistent scheme similar to 
the one of section 3.2J , i.e., using a variation principle for 
4>(R) and requiring that o~(R) balance the pressure differ- 
ence across the outer membrane. 

Discretization of the integral in equation (M) and tak- 
ing the variation of F with respect to <fi(R) give the bound- 
ary condition for the concentration gradient, 



F 



1 



dr-(j){x 2 [G(x)+xg(x)}cf> 2 + v} 



1 



+^KR i -Q(j)(R) 



<\4>{R) 



47T ~ 



2 rv ' dr 
and identify the integrand as the local pressure, 



P=\cf>{x 2 [G(x) 



xg(x)}4> 2 + /i}. 



(9) 



(10) 



Local balance between the pressure gradient and mem- 
brane tension is accounted for by a Laplace equation, 



1 



-Jiipi-i 



Pi 



1 nt \ d P 



(11) 



where &i is the rescaled tension (having dimension of length) 
of membrane i, ri its radius, and pi the pressure just out- 
side it. In the second equality we have assumed that the 
pressure profile is a smooth function on the length scale 
of D, the inter-membrane spacing. Equation (0) ensures 
that the pressure gradients resulting from the concentra- 
tion profile are consistent with the presumed tension pro- 
file. For the sake of mathematical convenience, we shall 
treat the tension profile as a smooth function as well, 
i.e., represent spatial variations of a by a first derivative, 
da/dr. Since the model allows for sharp changes in a, 
this approximation is merely technical and not physically 
corroborated. Consequently, the sharp changes will show 
up as singularities in an otherwise smooth tension profile, 
and will have to be treated separately to ensure that the 
smoothness of </> and p is maintained. When equation (nOf) 
is substituted in equation (11), the smoothness assump- 
tion leads to the following, second profile equation: 



[2G(x)+4xg(x)+x 2 g'(x)} 



-{3[G(x) + xg(x)} 



Ida; 
x dr 
M ,1 ' 



}- 



x 2 <p 2 (j> dr 



4f 
r 



0, (12) 



d(f>(R) 
dr 



a 



{x 2 [G{x)-xg{x)]4> 2 + 2x4> + n} 



r=R 



(13) 

Variation of F [Eq. (|9|)] with respect to R yields the me- 
chanical equilibrium (Laplace) equation for the outer mem- 
brane, 



_^ (i?) __ 



-R{p[^R),a(R)}-p }, (14) 



where p[4>(R), o~(R)] is the pressure just inside the outer 
sheet, whose dependence on (j)(R) and a(R) is given by 
equation (|lO|). The multiplier po, coupled to the total vol- 
ume [see equation (pp], is the external osmotic pressure 
exerted on the onion, i.e., the pressure just outside the 
outer sheet. Hence, the left-hand side of equation ( |l4|) is 
identified as the tension of the outer layer, a(R). This ob- 
servation, together with equation (O) and the definition 
of x, a = eSx 2 (f) 2 , lead to two equations for x(R) and <fi(R), 

x 2 [G{x) - xg(x)](j) 2 + 2(x + ex 2 )<f) + fi = (15) 



x 2 [G(x) + xg(x)]cj) 2 



AeS 



■x 4> + fi 



2po 



0. (16) 



One can proceed using boundary conditions (y_5|) and 
(|l6|). Yet, the formulation can be further simplified by 
employing an additional assumption. Prior to dilution, the 
onion is in contact with the outer membranes of neighbor- 
ing onions in the close-packed phase. This situation per- 
sists as long as the unbinding point of the regular L Q phase 
has not been reached. During this initial stage of dilution 
the inter-membrane spacing inside the onion is equal to 
the spacing between the outer membranes of neighboring 
onions (i.e., (j) is uniform throughout the sample). Conse- 
quently, the pressure is equal on both sides of the outer 
membrane and its tension vanishes, x(R) — 0. The profile 
equations (@) and (02) then have the trivial uniform solu- 
tion x(r) = 0, <j)(r) = 4>(R). When the unbinding point is 
reached, fi < /J.bin, the onion becomes separated from the 
surrounding membranes. The undulation pressure is still 
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exerted on the outer membrane from inside but vanishes 
outside, and tension must appear. Since the external pres- 
sure should become very low compared to the internal one, 
we may neglect po and assume that the internal pressure 
is balanced primarily by tension. This assumption leads 
to the following simplified boundary conditions: 



<fi(R) 



X (l + px 2 ) 



x 3 g(x) 



(17) 



r=R 



2[G{x) + xg(x)](l + px 2 



8eS 1 + px 2 



[x 2 g(x)]' 2 



XR xg(x) 



r=R 



_P_ 
Mbir 



(18) 



where p = (e/x)(l + 28/ R). 

Comparison between the boundary conditions (P-7|)- 
( |l8| ) and equation (||) shows that for jjl — > fj^ in the bound- 
ary values, x(R) and (p(R), coincide with the binodal ones, 
i.e., x(R) -> 0, <f>(R) -> 0bin, and dc/)(R)/dr -> 0. Thus, 
the outer layer continuously acquires different features as 
the membranes surrounding the onion unbind. The con- 
tinuous departure of x(R) from zero is also demonstrated 
in Fig. 0. The increase of x(R) when p. becomes slightly 
lower than p\,\ a is found from equation (Ffq) to scale like 



X(R) ~ \fX -Mbin| 



1/2 



M -» Mbi, 



(19) 



Hence, the tension in the outer layer, a(R) ~ x 2 (R), in- 
creases linearly with \p — /J.bm\- In addition, as p is further 
decreased, it is verified that x(R) of equation ( |l8| ) always 
remains smaller than the binodal value of equation (pi), 
required to stabilize bound membranes. The outer mem- 
branes of different onions, therefore, remain unbound from 
one another throughout the dilution. 



3.4 Profiles and Interfaces 

In order to calculate concentration and tension profiles 
in the onion, one should solve the profile equations (p) 



subject to the boundary conditions (13) 



and 

and flT8|). This is a system of coupled nonlinear equations, 
which is solved numerically. Note, however, that the en- 
tire formulation given above applies only in the case of low 
surfactant concentration (cf> <C 1). In order to correctly ac- 
count for concentrated domains in the onion (0 ~ 1) we 
actually solved a modified, more complicated set of equa- 
tions, as presented in the Appendix. Before considering 
the numerical results, it is useful to examine some general 
features of the solutions. 

Figure shows the x dependence of the local pres- 
sure for a given <f> [cf. Eq. (J1C|)] . There are two station- 
ary points of the pressure as a function of x: x = and 
x = x\ ~ 1.17. They correspond to two singularities of 
equation ( |l2| ) for the tension profile. Recall that the model 
allows for sharp changes in x(r), as long as the pressure 
profile remains smooth. The singularities in the tension 
profile signal such jumps as a consequence of the smooth- 
ing assumption in equation (|12|); they must be treated 



separately, as follows. If the profile reaches the singularity 
x(r) = x\, it cannot jump to any other value of x without 
violating the smoothness of p. Since equation ( fl2| ) is 1st- 
order, this implies that the rest of the profile must stay at 
X\. A different behavior is expected if the singular point 
x(r) = is reached. In this case the profile cannot remain, 
in general, at x = (zero tension), since this would imply 
also uniform pressure and, hence, a uniform concentration 
profile. The available option for the next membrane is to 
jump to a value of x = X2 — 1.91 while maintaining a 
continuous pressure profile. The meaning of the latter ob- 
servation is that a sharp interface in the tension profile, 
between a low-tension domain (x ~ 0) and a high-tension 
one (x ~ 1), is permitted. 

We now turn to the results of the numerical analysis. 
The analysis is restricted to stages of dilution where p 
has not reached very low values, so that x(R) could be 
assumed smaller than x\ (see, e.g., the values of x(R) in 
Fig. ||). The solutions to the profile equations are divided 
in this case into three families: 

(i) profiles which descend from x(R) to x — 0, then jump 
sharply to x — x-i and descend to remain at x = x\ 

(Fig- 1) 
(ii) profiles which, like family (i), have a sharp jump from 

x = to x = X2, but then proceed to higher values of 
z(Fig.g) 
(iii) smooth profiles which ascend from x(R) to x — X\ and 
remain there (Fig. [7|). 

Family (iii) may be called 'uniform', as it does not ex- 
hibit a sharp jump in tension. Indeed, since the pressure 
monotonously increases when going into the onion and x 
is uniform throughout the inner part of the onion, the 
concentration must increase as well and not remain uni- 
form. This increase, however, is only logarithmic in r [as 
can be found from equation (112) with constant x\. The 
'belt' family (i) is mostly uniform as well — it has only a 
narrow region of relatively high tension. Most interesting 
is the 'core' family (ii). In this case the jump in x di- 
vides the onion into an outer, low-tension region, and an 
inner, tense one. As we go inward, the tension continues 
to increase until diverging at a finite radius, where the 
membrane stack reaches the maximum concentration of 
close packing (4> = 1). (Note that the close-packing limit 
4> — > 1 does not imply a divergent free energy, since it is 
accompanied by a — > oo (see Fig. o). The high tension 
exponentially s uppr esses the appropriate term in the free 
energy [cf. Eq. (|AJ)].) 

The parameter space for numerical study is vast. Af- 
ter rescaling all distances with R, we are left with five pa- 
rameters: x, e, S, fl, and /j,. Nonetheless, the qualitative 
features, such as the three families of solutions described 
above, are found to be robust over a wide range of param- 
eter values. Figure g shows examples of 'profile diagrams' 
for three cuts through the parameter space. Dilution be- 
yond /Xbin usually results first in belt profiles. Cores may 
subsequently appear and, finally, the profile shifts to the 
uniform family for low enough fi. The structure of the di- 
agram, however, may be more complicated than this sim- 
ple sequence (Fig. ||B). Note the wide range of dilution 
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over which onion cores may be stable. Since fj, is not ex- 
pected in practice to become much lower than a few times 
Mbin, cores may remain stable over the entire process of 
dissolution (as indeed observed in experiments Jl^l). Not 
surprisingly, we find that core formation is promoted (i.e., 
occurs at higher /x) by lower membrane rigidity or higher 
temperature (smaller s ~ k/T), and stronger attraction 
between membranes (larger x)- 

The values assigned to Q in Fig. |] are rather large (a 
particularly large value was taken in Fig @B in order to 
demonstrate a richer diagram). Diagrams for smaller Q 
still exhibit the belt-to-core transition, yet spatial varia- 
tions occur on shorter distances from the boundary, re- 
sulting in bigger 'cores'. Another source of quantitative 
uncertainty is the factor b entering the definition of e 
[e = (n/T) / (wb)]. With b in the range between 0.06 (simu- 
lation Q) and 0.2 (theory [|o|), one needs k/T <~ 0.1- 
0.5 in order to get e <~ 0.5, as required for core forma- 
tion for reasonable values of \x (Fig. ||). These values are 
lower than the ones expected in the relevant experimental 
systems (« ~ a few T). It should be stressed, however, 
that we seek in this work merely qualitative mechanisms, 
rather than an accurate predictive capability. 

The narrow shell of high tension appearing in both the 
'belt' and 'core' profile families does not have, within the 
current approximation, a significant effect on the concen- 
tration profile (see Fig. H). Hence, unless an experiment 
is devised which will be sensitive to membrane tension, 
rather than density, this feature might not be directly 
observable in experiments. Nevertheless, the highly tense 
membranes may affect the kinetics of dissolution (e.g., 
hinder or assist water penetration and passage formation). 
This might explain the experimental observation of con- 
centric 'cracks' or 'rings' in diluted onions |15f|. 



an interface can be stabilized in a confined geometry with- 
out resorting to competing interactions. Stability becomes 
possible due to the tension in the membranes, which allows 
the system to avoid the usual equal-pressure condition of 
coexistence. In the absence of a length scale arising from 
competing interactions, the size of the inner domain must 
scale with the onion radius. The proportionality factor, 
however, may be small (see Fig. @). 

The conclusions drawn from the model rely on numer- 
ical solution of the profile equations and boundary condi- 
tions for specific values of parameters. Nevertheless, the 
general behavior presented in section ^J (e.g., the shift 
between the three families of profiles) is found to be ro- 
bust under change of parameters and even certain changes 
of the boundary conditions. Hence, we believe that the 
qualitative mechanisms indicated by this model are fairly 
general. The syste m of coupled nonlinear profile equations 
derived in section 3^ may, in principle, produce a much 
wider variety of solutions. It should be interesting, there- 
fore, to explore further (possibly less physical) areas of the 
parameter space. 

The theory presented here should be regarded as a pre- 
liminary step towards understanding the dissolution of the 
onion phase. In particular, it is focused on the first stages 
of dilution, where individual onions swell while maintain- 
ing their integrity. Further stages of dissolution involve 
onion breakage and coalescence, where intriguing instabil- 
ities are observed [13 . These phenomena probably require 



an altogether different theoretical approach. 
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4 Summary 

We have presented a theory for the swelling of the onion 
phase upon addition of solvent. A membrane stack in a 
spherical onion configuration can remain stable far be- 
yond the point where a flat lamellar phase disintegrates. 
This stability is achieved due to a tension profile acquired 
by the stack. As a result, the eventual dissolution of indi- 
vidual onions is bound to rely on membrane breakage and 
coalescence, thus taking very long time. 

At the unbinding point of the regular L Q phase, when 
individual onions become separated, a concentric shell of 
membranes with high tension ('belt') should first appear 
in the onion. If the membranes are not too rigid, subse- 
quent formation of a dense core might occur, as was ob- 
served in experiments. The cores may remain stable under 
extensive dilution (cf. Fig. g). In other cases, they may 
eventually disappear, giving rise to more uniform profiles. 

The interface formed between the inner, tense part of 
the onion, beyond the 'belt', and its outer part is ther- 
modynamically stable (on time scales shorter than hours) 
— the chemical potentials in the two domains are equal, 
and the pressure gradient is exactly balanced by an ap- 
propriate tension profile. This is a special example where 



Appendix: Expressions for High Concentration 



The various expressions derived in the previous sections 
apply only in the limit of low surfactant volume fraction, 
i.e., when the inter-membrane spacing is much larger than 
the membrane thickness, </> = S/D <C 1. Onion cores, how- 
ever, are regions of close-packed membranes, <j> ~ 1. Hence, 
a reliable calculation of profiles and profile diagrams, such 
as those presented in Figs. J5|-|3|, requires modified equa- 
tions, which are valid for high volume fractions as well. 
The following Appendix presents these modified expres- 
sions. 

The Helfrich interaction between tensionlcss membranes 
[ p0[ , Eq. (pi), is readily generalized to the case of finite 
membrane thickness. One uses the same arguments, yet 
the available space for undulations is now D — S instead 
of D. The resulting interaction energy per unit area is 



fund(D) 



bT 2 



n(D-8f 



(Al) 
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This leads to the following 'Flory-like' free energy density 
for a tensionless stack, 



/(</>) 



2(1 



x<t> - v4, 



(A2) 



which replaces equation (|2]). The modified expressions for 
the binodal and spinodal arising from this free energy are 
[compare to equation (pj)] 



Mbin = -(x/2)0bin(l + </>bin), 0bin(l - 4>bm) ' = X 

Ms P = -X0s P (l + 40 sp /3 - 03 P /3), 



sp (l-0 sp )- 4 = 2x/3. 



(A3) 



Seifert's calculation of the fluctuation-induced interac- 
tion in the presence of tension p7J is readily extended as 
well. The only modification required in equation (0) is the 
replacement of D with D — 5. The resulting 'Flory-like' 
free energy for a tense stack [generalizing Eq. (H)] is 



f(4>,a) = -x 2 G(x) 



(1 



x<t> - M 



(A4) 



with a modified definition of x, a; = (D—6)/It = (<V^t)(1- 



Writing a Ginzburg-Landau free energy similar to equa- 
tion (In) and taking the variation with respect to <f>{r < R), 
we obtain the modified version of the first profile equation 
[which replaces equation (g)], 

nV 2 cp+x 2 [G{x)-xg{x)/{l^<i ) )][cp/{l-cp)] 2 +2x<t>+^ = 0. 

(A5) 

Rewriting the free energy after minimization [cf. Eq. (H)], 

we identify the local pressure as 

p = ^{x 2 [G(x) + xg(x)/(l - 0)] [0/(1 - 0)] 2 + /x}, (A6) 

which replaces equation dlQ).n Substitu ting the modified 
expression for the local pressure, Eq. ( |A6|) , in the self- 
consistency condition, Eq. (|ll|) , we get the modified profile 
equation [compare to equation (|12|)], 

2(1 - (j>)G(x) + (4 - 4>)xg(x) + x 2 g'{x) dx 
x(l — (f>) dr 

[2 + (1 - 0) 3 ](1 - <f>)G{x) + [2 + (1 - 0) 2 ]x 5 (x) 



(1 - 0) 2 M 



2J.2 



x 4> 



(1-0)' 

a) dr r 



(A7) 



Finally, repeating the calculations described in sec- 
tion 3.3 , we find the following expressions for the boundary 



conditions 




G(x) 



xg(x) 



1- 



1- 



+ 2x0 + ft > 

(A8) 

1 Note that the stationary point of p as function of x for 
fixed (f>, defined in section pi as x = Xi (cf. Fig. W), is no longer 
a constant but depends on the value of </>. The same is true for 
x — X2, for which p(x — xz) — p(x = 0). 



which replaces equation (D_3j), and 

4> x[(i - <t>) 2 + px 2 ] 



(A9) 



1 — x 3 g(x) 

2[{l-4>)G{x)+xg{x)][(l-4>) 2 +px 2 ] 2 

[x 2 g(x)} 2 
_8^ (l-^+^ =JL 
X-K sg(a;) /Xbin 

instead of equations ( |l7| ) and (Jig) . 

The modified expressions derived in this Appendix were 
used to produce Figs. 
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Fig. 1. Formation of 'onion cores' upon dilution. (A) SDS-octanol-brine system J13| . (B) AOT-brine system [la. (Pictures 
courtesy of M. Buchanan.) 
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Fig. 2. Phase diagram of a lamellar phase with tension [Eq. (0)]. (A) Binodal and spinodal tension-chemical potential lines. 
(B) Values of x(fx) corresponding to the binodal. 
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Fig. 3. Boundary value of x as a function of chemical potential [Eq. (flq) with % — 0.1, £ = 0.1 and 8/R — 10 3 ]. 




Fig. 4. Dependence of pressure on x for given (f> [cf. Eq. (lid) 
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Fig. 5. 'Belt' profile exhibiting a narrow region of high tension. Solid line 

s/r = irr 3 , n/R 2 = o.oi, /u/Kml = -i.i. 



(r), dashed — x(r). Parameters: \ = 0.1, e = 0.1, 
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Fig. 6. 'Core' profile. Left — volume fraction (solid) and tension (dashed) profiles; note the logarithmic scale of position. Right 
— the volume fraction profile redrawn as a density plot. Parameters as in Fig. H except /J,/\fJlhin\ — —1-4 
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Fig. 7. 'Uniform' profile exhibiting a continuous tension profile. Parameters as in Fig. Ej except /ti/|/Ubin| = —6. 
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Fig. 8. 'Phase diagrams' of profile types encountered during onion swelling, (a) Dependence on e ~ n/T. Values of parameters 
are \ — 0.1, S/R — 10 -3 , and Q/R 2 — 0.01. (b) Same as (a) except Sl/R 2 = 0.1. (c) Dependence on \. Values of parameters 
are e = 0.21, S/R = 10~ 3 , and Q/R 2 = 0.01. 



